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MODELING  UNIDIRECTIONAL  COMPOSITE  LAMINATES  USING  XFEM 


INTRODUCTION 


Fiber-reinforced  polymers  (FRPs) — composites  consisting  of  a  thennoset  or 
thennoplastic  polymer  matrix  reinforced  by  carbon,  glass,  or  aramid  fibers — have  the  potential  to 
dramatically  influence  the  design  and  manufacturing  processes  of  innumerable  industries.  The 
significance  of  FRP  composites  is  that  they  can  be  customized  in  the  manufacturing  process  to 
meet  the  needs  of  specific  applications.  For  example,  FRPs  are  rapidly  becoming  the  preferred 
material  of  construction  for  applications  that  require  a  superior  strength-to-weight  ratio  and  a 
high  load-carrying  capacity.  FRPs  are  of  particular  interest  to  the  U.S.  Navy  because  of  their 
potential  use  in  underwater  applications  where  corrosion-resistance,  weight,  and 
implosion-resistance  are  indispensable  design  factors. 

Because  of  the  widespread  availability  of  measurement  techniques,  experimental  testing 
of  composite  materials  has  largely  outpaced  the  computational  modeling  ability  of  such 
complicated  materials — thus,  relegating  the  design  of  composite  structures  to  a  build-test-build 
methodology.  Recent  advancements  in  computational  physics-based  modeling  (PBM),  however, 
are  changing  the  manufacture  of  FRPs:  comprehensive  modeling  techniques  can  significantly 
reduce  the  effort  required  to  build  and  test  future  composite  structures  and  will  allow  previously 
untested  composite  arrangements  to  be  properly  analyzed  prior  to  any  physical  testing,  vastly 
reducing  cost  and  time  requirements. 

The  finite-element  method  (FEM)  has  been  a  popular  computational  analysis  technique 
for  modeling  mechanical  deformations  and  failure  for  any  material,  including  FRPs.  Numerical 
computation  of  FRP  fracture  has  historically  consisted  of  imposing  failure  conditions  on  a 
homogenous  material  with  averaged  composite  properties  to  capture  one  of  the  common  failure 
methods,  such  as  delamination.  Recent  research  has  investigated  fracture  using  the  extended 
finite-element  method  (XFEM).  Initiation  and  propagation  of  multiple  cracks  can  be  modeled 
simultaneously  without  mesh  refinement  or  successive  re-meshing  near  the  crack  tip.  Crack 
propagation  is  modeled  using  the  virtual  crack  closure  technique  (VCCT)  that  is  typically 
applied  to  linear  elastic  fracture  mechanics  (LEFM). 

This  report  documents  a  three-dimensional  (3-D)  computational  model  of  a  portion  of  a 
carbon  fiber/epoxy  composite  cylinder,  incorporating  both  cohesive  interface  failure  as  well  as 
the  XFEM  capability  available  in  Abaqus  v6.13.  A  traction-separation  cohesive  failure  model  is 
used  to  evaluate  failure  in  the  fiber/matrix  interface.  XFEM  is  used  to  model  crack  growth  and 
propagation  throughout  the  fibers  and  matrix.  The  cylinder  segment  is  loaded  under  external 
pressure  to  simulate  implosion  due  to  hydrostatic  depth  pressure.  Rather  than  modeling  the 
composite  as  a  single  material  (for  example,  homogenization),  individual  groups  of  fibers  are 
modeled  within  a  distinct  matrix  material.  The  composite  consists  of  either  one  layer  of 
unidirectional  (UD)  carbon  fibers  or  two  layers  with  0°/90°  fiber  orientation.  Three  groups  of 
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fibers  are  modeled  for  each  layer.  The  combination  of  both  cohesive  interface  failure  and 
material  failure  using  XFEM  in  one  model  allows  the  cohesive  failure  and  fracture  thresholds  to 
be  determined  for  given  inputs.  These  inputs  may  include  resin  and  fiber  material,  volume 
fraction,  fiber  layup,  interfacial  cohesive  strength,  and  material  fracture  toughness.  Analyses 
with  varying  input  parameters  are  presented  to  highlight  critical  thresholds  for  the  given 
geometry  and  loading  condition. 

With  the  infonnation  presented  in  this  report,  future  composite  structures  can  be  analyzed 
with  multiple  methods  of  failure  investigated  simultaneously.  Using  previously  obtained 
composite  material  properties,  thorough  modeling  of  future  designs  can  save  significant  time  and 
effort  in  experimentally  testing  and  modifying  composite  constructions. 


HISTORICAL  DEVELOPMENTS  IN  FRP  COMPOSITE  EVALUATION 


EXPERIMENTAL  TESTING 

FRP  composites  consist  of  a  thennoset  or  thennoplastic  polymer  matrix  reinforced  by 
carbon,  glass,  or  aramid  fibers.  While  individual  fibers  have  diameters  around  5  to  10  pm, 
thousands  of  fibers  are  typically  bundled  together  to  fonn  a  tow.  Single-direction-fiber-tow 
layers  can  be  (1)  stacked  into  varying  orientation  plies  or  (2)  they  can  be  woven  together.  The 
nature  of  manufacturing  FRPs  provides  a  unique  ability  to  customize  the  FRP  for  specific 
applications.  Fiber  and  matrix  material  selection,  fiber/matrix  volume  fraction,  and  laminate 
organization  are  just  a  few  of  the  many  factors  associated  with  designing  and  manufacturing  a 
desired  composite  material. 

Although  numerous  options  exist  when  the  use  and  design  of  FRP  composites  are  being 
contemplated,  the  exact  mechanical  nature  of  the  composite  has  not  been  fully  evaluated. 

Typical  research  of  composite  materials  has  consisted  of  extensive  testing  of  previously  built 
structures.  For  example,  composite  cylinders  designed  for  underwater  applications  are  first 
subjected  to  rigorous  impact  and  implosion  testing  to  detennine  their  structural  efficacy.1,2 
Tarpani  et  al.  determined  quasi-static  and  impact  properties  of  various  aeronautical-grade 
carbon/epoxy  laminate  structures  to  quantify  the  strength  of  various  composite  constructions. 

Furthennore,  the  effects  of  different  loading  scenarios  are  not  well  understood.  Although 
FRPs  show  particularly  good  material  properties  in  the  fiber  axial  direction,  transverse  loads  on 
complex  structures  can  lead  to  failure  at  lower  stress  states.  The  failure  of  FRPs  under 
compression,  which  leads  to  the  creation  of  fiber  kink  bands,  was  experimentally  studied  by 
Kyriakides  et  al.4.  Other  factors,  such  as  sensitivity  to  rate  of  defonnation,  can  be  measured  for 
various  composite  combinations.  ’  Kepple  et  al.  showed  that  novel  composite  organizations, 
such  as  carbon  nanotubes  grown  onto  woven  carbon  fiber  composites,  can  be  created  to  combat 
the  shortfalls  of  traditional  composite  structures. 
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FRP  FAILURE  METHODS 


Much  work  has  been  done  to  attempt  to  comprehensively  analyze  and  model  the 
constituents  of  a  composite  material.  FRPs  fail  by  one  or  more  of  four  different  methods:  (1) 
fiber  fracture,  (2)  matrix  cracking,  (3)  delamination,  or  (4)  fiber  pull-out.  The  most  common 
causes  of  these  failure  methods  are  the  failure  of  the  matrix  (matrix  cracking,  delamination)  or 
the  failure  of  the  fiber/matrix  interface  (fiber  pull-out,  delamination).  Each  failure  phenomenon 
corresponds  to  a  failure  in  one  of  three  primary  modes  of  material  fracture:  crack  opening, 
in-plane  shearing,  or  out-of-plane  tearing.  A  review  of  the  various  modes  of  fracture,  which  are 
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discussed  in  the  next  section,  is  provided  by  Prasad  et  al.  A  comprehensive  overview  of  the 
mathematical  modeling  of  different  fracture  phenomena  in  fiber  composites  is  provided  by 
Mishnaevsky  and  Brondsted.9  In  addition,  Pinho  et  al.10  developed  physically  based  criteria  for 
calculating  multiple  modes  of  failure  in  FRPs. 

These  composite  structures  consist  of  complex  material  interactions  and  can  be  analyzed 
at  differently  sized  scales.  The  majority  of  research  has  focused  on  investigating  composites  in  a 
large  scale,  where  the  composite  is  assumed  to  be  a  homogenous  material  with  averaged 
mechanical  properties.  Other  research  has  focused  on  the  random  nature  of  individual  fibers  in 
the  microscopic  scale.9’ 11 
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Fracture  of  composite  laminates  has  been  modeled  using  a  variety  of  FEMs.  Krueger 

documented  a  large  review  of  fracture  toughness  characterization  for  multiple  modes  of  fracture 

using  the  VCCT.  Sun  et  al.  ~  developed  an  FEM  using  plate  elements  to  characterize  fracture 

toughness  in  woven  carbon/epoxy  composites  under  mixed-mode  loading.  Waseem  and 

Kumar14  investigated  delamination  of  carbon  fiber/epoxy  composites  using  multiple  fracture 

techniques.  Camanho  et  al.15  used  interfacial  decohesion  elements  between  layers  to  model 

delamination.  Delamination  of  materials  has  also  been  studied  by  Zhao,16  who  focused  on  the 

viscoelastic  behavior  of  composites  under  elevated  temperature.  Wiinmer  et  al.  analyzed  onset 

and  growth  of  cracks  in  laminated  composites  in  curved  structures.  With  rigorous  modeling 

efforts,  certain  capabilities  of  existing  material  constructions  can  be  improved  upon.  Segala  and 
18.. 

Cavallaro  investigated  the  failure  modes  of  UD  fibrous  composites  subjected  to  ballistic 
impacts  in  an  effort  to  increase  energy  absorption  capacities. 

While  most  of  the  research  on  the  fracture  of  composite  materials  has  used  conventional 
FEM  techniques,  the  XFEM  provides  unique  capabilities  in  modeling  fracture.  Though  standard 
finite-element  analyses  require  a  highly  refined  mesh  near  the  assumed  crack  location,  XFEM 
allows  cracks  to  fonn  and  propagate  along  a  natural  crack  path  without  prior  knowledge  of  its 
location.  In  addition,  continual  re-meshing  near  the  crack  front  is  not  required.  Eliminating 
these  computationally  expensive  methods  makes  modeling  fracture  using  XFEM  markedly 
easier.  Widespread  research  has  been  conducted  implementing  XFEM  into  established  finite- 
element  codes.  Recently,  commercial  software  packages  such  as  Abaqus  have  included  an 
XFEM  capability.  With  these  new  abilities,  previously  impractical  simulations  involving 
composite  structures  can  be  analyzed  quickly  and  efficiently. 
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Primarily,  XFEM  can  be  used  to  validate  experimental  procedures  and  predict  future  test 
results.  Motamedi  and  Milani19  developed  a  3-D  model  utilizing  XFEM  and  the  cohesive  zone 
method  (CZM)  to  analyze  delamination  under  mode  I  loading  (modes  I,  II,  and  III  are  discussed 
in  the  Fracture  Mechanics  section).  Leven  and  Rickert20  investigated  the  integrity  of  modeling 
various  cracks  using  XFEM.  Grogan  et  al.  used  an  XFEM-based  methodology  to  validate 
mode  I  and  II  fracture  when  analyzing  thennal  fatigue  of  carbon-fiber-reinforced  polymers 
(CFRP). 

With  the  XFEM,  a  wide  range  of  composite  materials  can  be  evaluated.  For  example, 
Benvenuti  et  al.22  applied  XFEM  concepts  to  modeling  FRP -reinforced  concrete  and  Sosa  and 
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Karapurath"  used  XFEM  to  model  delamination  of  fiber-metal  laminate  composites  with 
orthotropic  material  properties.  Moreno  et  al.24  modeled  and  analyzed  crack  propagation  in 
chopped  glass-fiber  composites  under  biaxial  loading  and  compared  the  results  to  experimental 

.95 

findings.  Motamedi  investigated  the  nonuniform  nature  of  composites  using  a  combination  of 
XFEM  and  other  commercial  FEM  software.  Instead  of  using  typical  isotropic  material 
properties,  Nagashima  and  Suemasu26  developed  an  orthotropic  material  model  undergoing 
mode  I  delamination. 


MODELING  INTERFACE  INTERACTIONS  BETWEEN  DISTINCT 
FIBER  GROUPS  AND  THE  MATRIX 


Although  extensive  research  has  been  conducted  in  modeling  FRPs  using  finite-element 
techniques  at  both  micro  (individual  fibers)  and  macro  (homogenous  composite  materials)  scales, 
little  work  has  been  done  to  model  interface  interactions  between  distinct  fiber  groups  and  the 
matrix  in  conjunction  with  matrix  fracture  using  XFEM.  A  goal  of  this  research  effort  was  to 
develop  a  three-dimensional  (3-D)  finite-element  model  that  incorporated  multiple  failure 
techniques  available  in  Abaqus  v6.13.  This  research  focused  on  separating  individual  tows  of 
fibers  from  the  surrounding  matrix.  Distinct  cylindrical  groups  of  fibers  were  treated  as  a 
homogenous  material  with  a  smooth  interface  with  the  surrounding  matrix.  By  defining  the 
model  in  this  way,  interfacial  cohesion  could  be  implemented  while  the  overall  cylindrical 
geometry  was  maintained.  In  addition  to  evaluating  cohesive  interface  failure  between  the  fibers 
and  matrix,  this  research  also  evaluated  matrix/fiber  fracture  using  the  XFEM. 

The  developed  model  was  used  to  study  the  onset  of  cracking  and  decohesion  in 
unidirectional  one-  and  two-ply  carbon-fiber/epoxy  composite  cylinders.  The  loading  condition 
on  the  composite  cylinder  was  exterior  hydrostatic  pressure — a  particularly  useful  loading 
condition  because  of  its  widespread  effect  on  a  variety  of  underwater  systems.  Compression 
caused  by  external  pressure  results  in  material  fracture  and  eventual  implosion,  presenting  a 
significant  concern  for  designers. 
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FRACTURE  MECHANICS 


In  a  generic  material,  three  different  modes  of  fracture  can  occur  (see  figure  1): 

1 .  Mode  I  fracture  results  in  two  crack  faces  being  pulled  away  from  each  other. 

2.  Mode  II  fracture  occurs  because  of  in-plane  shear  loading,  where  the  two  fractured 
surfaces  slide  with  respect  to  each  other  along  the  crack  plane. 

3.  Mode  III  fracture,  out-of-plane  shear  due  to  shear  loading,  occurs  when  the  crack 
faces  separate  out  of  the  crack  plane. 

The  three  failure  modes  illustrated  in  figure  1  correspond  to  the  various  modes  of  failure 
in  FRPs  previously  mentioned:  fiber  fracture,  matrix  cracking,  delamination,  or  fiber  pull-out. 
Delamination  is  the  result  of  a  mode  I  fracture  where  separate  layers  of  a  composite  are  tom 
away.  Fiber  fracture  and  matrix  cracking  are  typically  a  function  of  mode  II  or  III  fractures 
where  the  material  undergoes  shear  failure.  Fiber  pull-out  happens  when  the  cohesive  interface 
between  the  fiber  and  matrix  severs,  allowing  the  fiber  to  slide  through  the  matrix. 


(a)  (b)  (c) 


Opening  In-Plane  Shear  Out-of-Plane  Shear 

Figure  1.  (a)  Mode  I,  (b)  Mode  II,  (c)  Mode  III  Fracture 


The  onset  and  growth  of  failure  can  be  determined  by  comparing  the  strain  energy  release 
rate  Gy  to  a  critical  fracture  toughness  value  Gc ■  Gc  is  an  experimentally  determined  material 
property  and  can  be  broken  down  into  a  combination  of  the  individual  mode  fracture  toughness 
values  Gjc,  Gnc,  and  GUic-  GT  is  a  summation  of  Gj,  Gn,  and  Gin,  which  correspond  to  energy 
release  rates  due  to  mode  I,  II,  and  III  fractures,  respectively.  When  Gt  >  Gc,  the  crack  can  be 
expected  to  propagate.  Using  these  individual  strain  energy  release  rates,  an  overall  failure 
criterion  can  be  established.  Benzeggagh  and  Kenane  created  a  simple  mathematical  relationship 
combining  the  various  mode  strain  energy  release  rates  to  detennine  Gc,  which  is  shown  in 
equation  (1): 

-  GiC  +  (G//c  +  Gmc  -  GIC )  (  " Gt 

GT  -  G,  +  G, /  +  GUI.  (1) 
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Typically,  mode  I  fracture  toughness  Gic  is  determined  through  the  double  cantilever 
beam  (DCB)  test.  ’  ~  One  end  of  a  cantilever  beam  is  pulled  in  opposite  directions,  splitting  the 
beam  in  two.  One  of  the  most  common  experimental  methods  to  measure  mode  II  fracture 
toughness  Guc  is  the  end-notched-flexure  (ENF)  test.  The  test  places  a  pre-cracked  specimen 
into  a  three-point  bending  fixture.  The  crack,  located  at  the  end  of  the  specimen,  propagates 
along  the  length  of  the  beam,  resulting  in  delamination.  A  number  of  experiments  have  been 
proposed  to  measure  mode  III  fracture  energy  release  rate  Gmc •  One  such  method,  described  by 
Lee,29  used  the  edge-crack-torsion  (EOT)  technique  to  investigate  edge  delaminations. 

The  vast  majority  of  loading  conditions  on  composite  materials,  however,  involve 
mixed-mode  loading.  The  mixed-mode-bending  (MMB)  test  has  been  designed  to  encapsulate 
both  modes  I  and  II  fracture  simultaneously.  The  distance  of  the  lever  ann  can  be  modified  to 
provide  loading  that  corresponds  to  a  spectrum  between  pure  mode  I  loading  to  pure  mode  II 
loading.  An  example  of  the  mixed-mode  bending  is  given  by  Prasad  et  al. 


XFEM 


30 

The  XFEM  was  originally  developed  by  Belytschko  and  Black  and  is  based  on  the 
partition-of-unity  method  (PUM)  developed  by  Melenk  and  Babuska.31  In  PUM,  additional 
enrichment  functions  are  inserted  into  the  conventional  shape  functions  to  account  for 
discontinuities  across  the  crack  front.  This  continuity  allows  cracks  to  form  and  propagate 
without  prior  knowledge  of  the  crack  location  and  does  not  require  remeshing  throughout  the 
simulation  to  maintain  a  fine  mesh  near  the  crack  tip.  The  partition  of  unity  is  further  described 
in  detail  by  Belytschko  et  al.32  The  applied  enrichment  functions  introduce  additional  degrees  of 
freedom  to  affected  elements,  allowing  the  element  to  undergo  transformation  while  maintaining 
continuity.  In  fracture  mechanics,  the  enrichment  functions  consist  of  a  near-tip  asymptotic 
function,  which  captures  the  stress  singularity  at  the  crack  tip,  and  a  discontinuous  function  to 
map  the  displacement  of  the  crack  faces  within  an  element.  The  enriched  shape  function  is 
given  in  equation  (2): 

uh(x)  =  Y.11  N,  (x)u,  +  N,  (x)  [«(/(*))  -  H  (/(*;))]  q° 

+  Ij Zkzsc  NkO)  [VA*)  -  (2) 

The  first  summation  includes  the  standard  shape  functions  and  nodal  degrees  of  freedom, 
where  Nf  are  the  standard  shape  functions  and  uj  are  nodal  degrees  of  freedom.  The  second 
summation  adds  additional  enrichments  to  the  elements  in  the  set  Sh  encompassing  the  crack  tip. 
H  is  the  Heaviside  function  and  q j  are  enrichment  coefficients.  The  third  summation  adds 
enrichments  to  the  elements  cut  by  the  crack,  excluding  the  elements  inherent  in  the  second  tenn 
and  are  included  in  the  set  Sc-  The  tenn  is  a  set  of  enrichment  functions  near  the  crack  tip. 
H  is  defined  in  equation  (3)  as: 

=  {o!  otherwise  ‘  (3) 
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Because  the  crack  affects  only  elements  through  which  it  passes,  only  the  connecting 
nodes  gain  the  additional  degrees  of  freedom.  Figure  2  shows  a  schematic  of  how  the  nodes  are 
assigned  around  a  given  crack. 
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Figure  2.  Crack  Running  Through  a  Set  of  Elements  with  Enriched  Functions 


Material  Fracture 

In  Abaqus,  two  predominant  methods  can  be  used  to  calculate  the  movement  of  a  crack 
through  a  material:  one  method  is  based  on  the  cohesive  traction-separation  interaction,  and  the 
other  is  based  on  LEFM.  Because  the  purpose  of  this  research  was  to  investigate  fracture  of  a 
brittle  composite  matrix,  LEFM  was  the  chosen  failure  method.  The  stress  energy  release  rates 
were  determined  using  the  VCCT.  In  VCCT,  the  energy  absorbed  through  material  fracture  is 
assumed  to  be  the  work  required  to  close  the  crack  faces.  Crack  initiation  occurs  when  the 
prescribed  maximum  principal  stress  ops  is  reached.  Upon  initiation  and  growth,  the  crack 
follows  the  nonnal  of  the  maximum  tangential  stress  direction.  The  fracture  energy  release  rate 
is  calculated  using  the  Benzeggagh-Kenane  (BK)  law.'  The  mode  fracture  toughness  values  can 
be  calculated  experimentally  and  used  as  inputs  in  the  finite-element  model.  Table  1  lists  a  set  of 

25 

material  properties  typically  seen  in  fiber/epoxy  composites. 


Table  1.  Material  Properties  for  Common  Fibers  and  Epoxy  Matrix25 


Property 

E  msi  (GPa) 

V 

crv  ksi  (MPa) 

Carbon  Fiber 

15-60  (103-414) 

0.33 

250-350(1723-2413) 

Glass  Fiber 

10-12  (69-83) 

0.33 

500-650  (3450-4480) 

Kevlar  (Aramid) 

12-19(83-131) 

0.3 

406 (2800) 

Epoxy  Matrix 

0.35-0.7  (2. 4-4. 8) 

0.28 

7-15 (48-103) 
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COHESIVE  INTERFACE 


In  addition  to  matrix  fracture,  this  research  also  investigated  the  separation  of  the 
fiber/matrix  interface,  another  method  of  failure.  This  separation  occurs  when  shear  or  normal 
stresses  overcome  the  interfacial  bond,  resulting  in  the  fiber  sliding  through  the  matrix  surface. 

In  Abaqus,  cohesion  can  be  modeled  between  two  surfaces  in  contact.  Cohesive  surfaces  are 
assumed  to  be  thin  films  with  zero  thickness  that  fonn  a  rigid  bond.  To  break  the  cohesive  bond, 
a  damage  criterion  must  be  met.  Prior  to  damage  initiation,  the  two  surfaces  (fiber  and  matrix) 
possess  a  cohesive  stiffness  in  each  direction.  These  stiffness  values,  K„,  Ks,  and  Kh  represent  a 
linear  relationship  between  the  cohesive  stress  tn(ts,  tt )  and  cohesive  separation  8n(8s,  8t)  for 
the  normal  and  shear  directions.  The  damage  criterion  used  in  this  model  uses  the  maximum 
nominal  stress  in  either  the  nonnal  or  shear  directions,  shown  in  equation  (4): 


where  t°,  t°,  or  are  the  material  input  maximum  nominal  stresses  in  the  normal  and  shear 
directions,  respectively.  Upon  reaching  the  prescribed  maximum  stress  ),  it  is  assumed 

that  damage  has  occurred  in  the  cohesive  bond.  After  damage  initiation  has  occurred,  the 

cohesive  interface  follows  a  linear  damage  law  until  reaching  the  final  displacement  8 1  (d/,  8[). 
Figure  3  shows  the  damage  initiation  and  evolution  path. 


t 


The  overall  damage  in  the  linear  damage  law  is  defined  by  the  scalar  variable  D.  When 
D  reaches  a  value  of  1 ,  the  bond  is  assumed  to  be  completely  severed.  The  nominal  stress 
components  are  associated  with  D  as: 

t  =  f(!  -  D)Q,  tn>  0 
n  l  tn>  otherwise 

ts  =  (1  -  D)ts 

tt  =  (1  -  D)tt.  (5) 
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3-D  FINITE-ELEMENT  MODEL 


DESCRIPTION 

In  this  analysis,  a  cylindrical  geometry  was  chosen  because  of  its  widespread  use  as  an 
underwater  vessel.  Because  of  the  symmetrical  nature  of  cylinders,  only  a  slice  of  the  cylinder 
was  required  to  be  modeled;  therefore,  l/20th  of  the  cylinder  circumference  was  used  as  a  unit 
cell  of  the  entire  cylinder  within  a  cylindrical  coordinate  system.  Symmetry  boundary  conditions 
were  imposed  on  the  0-  and  axial  z-directions  to  both  the  fiber  and  matrix.  With  these 
dimensions,  three  fiber  groups,  or  tows,  with  a  diameter  of  0.02  inch  (0.5  mm)  were  modeled  for 
each  layer,  with  a  0.02-inch  separation  between  the  fibers.  The  two  layers  were  oriented  in  the 
0°/90°  directions,  with  0°  corresponding  to  the  cylinder  axis  and  90°  in  the  circumferential 
0-direction. 

This  orientation  was  chosen  to  highlight  the  differences  in  stress  and  fracture  in  the 
cylinder  warp  and  weft  directions.  The  cylinder  wall  thickness  was  0.08  inch  (2  mm),  allowing 
0.01  inch  (0.25  mm)  of  matrix  between  the  top  and  bottom  plies.  The  selected  geometry  was 
chosen  to  accentuate  the  matrix  filling  between  individual  carbon  fiber  groups,  allowing  cracks 
to  form  and  propagate  throughout  the  matrix.  In  addition  to  having  two  layers  of  fibers, 
single-ply  composites  with  fibers  in  the  axial  z-direction  were  analyzed  separately.  The  purpose 
of  analyzing  single-ply  composites  was  to  isolate  the  potential  damage  mechanisms  in  the  top 
layer,  as  damage  was  expected  to  occur  in  the  matrix  near  the  bottom  0-direction  fibers  first. 
Figure  4  shows  an  initial  configuration  of  the  two-ply  and  single-ply  models. 


Figure  4.  Two-Ply  (Left)  and  One-Ply  (Right)  Finite-Element  Model 


The  fibers  were  modeled  as  a  linear  elastic  orthotropic  material.  Damage  was 
incorporated  into  the  fibers  with  a  maximum  principal  stress  condition.  Above  the  maximum 
principal  stress  ops,  the  fibers  undergo  fracture  according  to  the  described  failure  model.  The 
matrix  was  modeled  as  an  isotropic  linear  elastic  material,  with  a  similarly  defined  maximum 
principal  stress  failure  criterion.  Tables  2  and  3  show  the  material  properties  used  for  both  the 
fiber  and  matrix. 
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Table  2.  Constant  Material  Properties  for  the  Fiber  and  Matrix 


Property 

En 

ksi 

(GPa) 

E22 

ksi 

(GPa) 

E33 

ksi 

(GPa) 

G12 

ksi 

(GPa) 

C,3 

ksi 

(GPa) 

G23 

ksi 

(GPa) 

V12 

Vn 

V23 

Fiber 

18500 

(128) 

1160 

(8) 

1160 

(8) 

435 

(3) 

435 

(3) 

667 

(4.6) 

0.28 

0.28 

0.46 

Matrix 

350 

(2.4) 

- 

- 

- 

- 

- 

0.3 

- 

- 

Table  3.  Fracture  Toughness  Values  for  the  Fiber  and  Matrix 


Property 

Gic 

lb/in  (N/m) 

Gnc 

lb/in  (N/m) 

Guic 

lb/in  (N/m) 

Fiber 

5 

5 

5 

(875) 

(875) 

(875) 

Matrix 

5 

5 

5 

(875) 

(875) 

(875) 

Using  XFEM  enabled  the  modeling  of  fracture  for  both  the  libers  and  the  matrix. 
However,  because  the  fibers  were  significantly  stronger  than  the  matrix,  and  because  of  the 
compressive  stresses  present,  only  the  matrix  was  expected  to  fail.  The  interface  between  the 
matrix  and  fibers  was  modeled  as  a  cohesive  surface  with  a  traction-separation  damage  law.  The 
cohesion  fails  when  the  nonnal  or  shear  stress  reaches  a  maximum  nominal  stress  tf ). 

After  damage  initiation,  the  traction-separation  law  activates,  giving  a  progressive  damage  that 

f  f  f 

reaches  full  value  at  the  specified  displacement  ( 8JS  ,8t).  Once  a  bond  initially  breaks,  it  will 

not  reconnect  with  subsequent  contact. 

In  order  to  investigate  the  effects  of  combining  element  fracture  and  cohesive  failure  into 
one  model,  the  values  for  t°)  were  varied  (see  table  4).  Cohesive  strength  is  evaluated  at 
t°(t° ,  t°)  =  1,5,  50,  and  500  ksi.  Maximum  principal  stress  in  the  matrix  required  for  crack 
forming  was  investigated  at  er^  =  0.5,  1,  5,  and  10  ksi.  A  static  analysis  was  conducted  with 
external  pressure  as  the  loading  condition. 


Table  4.  Fiber,  Matrix,  and  Fiber/Matrix  Cohesive  Stress  at  Failure 


f 

Fiber  a’vs 

Matrix  o™s 

Cohesive  t°,  tQs, 

ksi  (MPa) 

ksi  (MPa) 

ksi  (MPa) 

250 

0.5,  1,5,  10 

1,5,50,  500 

(1723) 

(3.45,  6.9,  34.5,  69) 

(6.9,  34.5,  345,  3450) 
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RESULTS 


Initially,  the  composite  cylinder  was  simulated  without  fracture  in  the  matrix  or  fibers 
and  with  a  sufficiently  high  cohesive  strength  to  prevent  any  cohesive  damage  from 
accumulating.  Figure  5  shows  the  Von  Mises  and  maximum  principal  stresses  in  the  matrix 
alone.  The  highest  Von  Mises  stress  is  seen  between  the  axial  fibers.  As  the  cylinder  is 
compressed,  large  compressive  stresses  are  present  in  the  circumferential  (9-direction.  With  the 
fibers  in  the  top  ply  aligned  perpendicular  to  the  highest  stress  direction,  the  surrounding  matrix 
must  support  the  load.  As  for  the  maximum  principal  stress,  the  inner  section  of  the  cylinder 
experiences  the  highest  stress  state,  beneath  the  circumferential  fibers.  It  is  at  this  location  that 
matrix  cracking  can  be  expected  to  form. 


S,  Mises 
(Avg:  759b) 

B+1.300e+01 
+  1.242e+01 
+  1.183e+01 
+  1.125e+01 
+  1.067e+01 
+  1.008e+01 
+  9.500e+00 
+  8.91 7e+00 
+  8.333e+00 
+  7.750e+00 
+  7.167e+00 
+  6.583e+00 
+  6.000e+00 


r 


S,  Max.  Principal 
(Avg:  759b) 

B+l.OOOe+OO 
+  4.167e01 
-1.667e-01 
-7.500e-01 
-1.333e+00 
-1.917e+00 
-2.500e+00 
-3.083e+00 
-3.667e+00 
-4.250e+00 
-4.833e+00 
-5.41 7e+00 
-6.000e+00 


(a) 


(b) 


Figure  5.  Stresses  in  Matrix  Without  Failure:  (a)  Von  Mises  and  (b)  Maximum  Principal 
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The  matrix  is  a  candidate  for  future  crack  growth;  the  fibers  are  responsible  for  carrying 
the  compressive  load.  Figure  6  shows  the  Von  Mises  and  maximum  principal  stresses  in  the 
fibers.  A  significantly  higher  Von  Mises  stress  state  (figure  6(a))  can  be  seen  in  the 
circumferential  fibers,  which  is  attributed  to  the  fact  that  the  primary  deformation  is  the  decrease 
of  the  cylinder  diameter.  Though  the  axial  fibers  experience  transverse  compression,  little  load 
is  supported  in  the  fiber  direction.  It  can  be  seen  in  figure  6(b)  that  the  fibers  in  a  maximum 
principal  stress  condition  are,  as  expected,  entirely  in  compression.  Fracture,  only  occurring  in 
tension,  is  not  expected  to  be  seen  in  the  fibers. 


S,  Mises 
(Avg:  754t>) 

+  5.750e+02 
+  5.279e+02 
+  4.808e  +  02 
+  4.338e+02 
+  3.867e+02 
+  3.396e  +  02 
+  2.925e  +  02 
+  2.454e  +  02 
+ 1 .983e  +  02 
+  1.5l3e  +  02 
+  1.042e  +  02 
+  5.708e  +  01 
+  1.000e+0l 


S,  Max.  Principal 
(Avg:  75^i) ) 

B+0.000e+00 
-5.000e-01 
-l.OOOe+OO 
-1.500e+00 
-2.000e+00 
■2.500e+00 
-3.000e+00 
-3.500e+00 
-4.000e+00 
-4.500e+00 
-5.000e+00 
-5.500e+00 
-6.000e+00 


(a) 


(b) 


Figure  6.  Stresses  in  Fibers  Without  Failure:  (a)  Von  Mises  and  (b)  Maximum  Principal 
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In  order  for  cracks  to  form  and  propagate  in  the  matrix,  the  maximum  principal  stress 
must  overcome  cr^.  In  the  interest  of  consistency,  rather  than  change  the  external  pressure  pe, 
a£l  is  reduced  to  induce  failure.  Figure  7  shows  the  translucent  final  configurations  with 
aJU  =  0.5,  1,  5,  and  10  ksi,  respectively.  When  fracture  occurs,  denoted  in  red,  it  is  located  near 
the  bottom  portion  of  the  center  circumferential  fiber.  Fracture  initiates  near  multiple  sites  along 
the  fiber,  then  propagates  in  the  circumferential  direction  until  reaching  the  boundaries. 
Presumably,  because  of  the  cylindrical  boundary  conditions  applied  to  the  model,  fracture  can 
occur  near  any  fiber  in  the  bottom  layer.  It  is  assumed  that,  after  initial  damage  accumulation, 
the  stress  state  is  relaxed  to  a  point  where  subsequent  fracture  sites  are  unlikely  to  fonn.  At 
o™s  =  0.5  ksi,  the  crack  propagates  along  the  entire  length  of  the  fiber.  The  crack  initiates  in 
multiple  sites  when  a£l  =  1  ksi  but  does  not  propagate  fully.  At  er™  =  5,10  ksi,  the  requisite 
failure  criterion  is  not  reached  and  no  cracks  are  fonned. 


Figure  7.  Crack  Formation  and  Propagation  in  ops  =  0.5  ksi  (I),  1  ksi  (II), 

5  ksi  (III),  and  10  ksi  (IV) 
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To  investigate  the  cohesive  interface  failure,  cr£l  was  maintained  at  5  ksi,  higher  than  the 
accumulated  principal  stress  in  the  model.  The  required  nominal  stress  to  initiate  cohesive 
failure  was  varied  with  t°)  =  1,  5,  50  and  500  ksi.  Figure  8  shows  the  cohesive  damage  D 
for  each  of  these  analyses.  At  t  °  =  1  ksi,  nearly  the  entire  interface  is  severed.  The  majority  of 
the  interface  was  damaged  with  t  °  =  5  ksi,  while  damage  occurred  only  in  a  few  locations  when 
=  50  ksi.  As  expected,  no  damage  occurred  when  =  500  ksi.  When  only  partial  damage 
was  present,  the  cohesive  bonds  initially  failed  at  the  diagonal  regions  of  the  cylindrical  fibers. 


Ill 


Figure  8.  Cohesive  Damage  D  with  (t°,  )  =  1  ksi  (I),  5  ksi  (II), 

50  ksi  (III),  and  500  ksi  (IV)  ksi 
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When  cohesive  interfacial  damage  was  combined  with  matrix  cracking,  one  can  see 
multiple  methods  of  material  fracture  in  a  single  model.  In  order  to  properly  assess  the 
combination  of  these  two  failure  phenomena,  the  matrix  maximum  principal  stress  prior  to 
failure  was  cr^  =  1  ksi  and  the  cohesive  interface  maximum  nominal  stress  was  t %  =  5  ksi. 

Figure  9  illustrates  the  Von  Mises  stress,  interfacial  damage  D,  and  the  crack  location.  Similar  to 
the  damage  shown  in  figure  8,  the  fiber/matrix  interface  is  partially  broken  around  all  fibers.  A 
crack  emerges  along  the  center  circumferential  fiber  and  grows  in  the  ^-direction. 
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Figure  9.  Cohesive  Damage  D  and  Crack  Formation  with  tn(t°,  t°)  =  5  ksi  and  a™s  =  1  ksi 
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In  the  previous  simulations,  cracks  formed  and  grew  along  the  circumferential  0-direction 
fibers;  moreover,  these  fibers  experienced  the  highest  stress.  If  the  0-direction  fibers  were 
removed,  the  structural  capacity  would  be  supported  by  the  matrix  and  transverse  compression  of 
the  z-direction  fibers.  When  the  circumferential  fibers  were  removed,  with  just  one  ply  of  axial 
fibers  in  the  cylinder,  fracture  occurred  in  the  matrix  along  the  direction  of  the  fibers.  This 
fracture,  however,  occurred  at  a  significantly  lower  external  pressure  pe  =  1  ksi.  Without  the 
reinforcing  circumferential  fibers,  the  matrix  undergoes  more  deformation  and  eventual  failure. 
Furthermore,  some  crack  propagation  was  observed  in  the  matrix  in  the  0-direction.  Along  with 
the  fonnation  of  a  crack,  the  cohesive  interface  was  more  severely  damaged.  Figure  10  shows 
the  failure  at  pe  =  1  ksi  for  a  representative  one-ply  simulation. 


Figure  10.  Cohesive  Damage  D  and  Crack  Formation 
with  £j)  =  50  ksi  and  o™s  =0.5  ksi 
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CONCLUSION 


In  this  report,  the  application  of  the  extended  finite-element  method  (XFEM)  to  a  3-D 
finite-element  model  was  combined  with  cohesive  interface  failure  mechanisms  to  capture 
multiple  methods  of  fracture  in  carbon  fiber/epoxy  composite  cylinders.  External  pressure  was 
used  as  a  loading  condition  because  typical  applications  of  composite  cylinders  as  underwater 
vessels  are  exposed  to  hydrostatic  depth  pressure.  In  the  interest  of  consistency,  geometry  and 
pressure  amplitude  were  held  constant  while  material  properties  were  modified.  Thresholds  for 
material  failure  prior  to  crack  initiation  were  varied  to  detennine  failure  thresholds  for  the  given 
geometry  and  pressure  loading.  In  order  for  fracture  to  occur  in  the  matrix,  the  maximum 
principal  stress  was  required  to  surpass  the  material  property  er^.  Fracture  was  not  seen  in  the 
fibers  because  that  stress  condition  was  entirely  compressive.  At  the  lowest  value  of  er^,  the 
crack  initially  formed  near  the  circumferential  ^-direction  fibers,  followed  by  propagation  along 
the  fiber  direction.  At  higher  a™s,  only  partial  cracks  fonned,  if  at  all. 

In  addition  to  modeling  crack  formation  and  growth,  cohesive  interface  failure  was 
investigated  for  varying  cohesive  strengths  £°.  Complete  cohesive  failure  occurred  at  the  lowest 
,  and  as  increased,  progressively  less  damage  accumulated.  When  only  partial  damage  was 
present,  the  cohesive  bonds  initially  failed  at  the  diagonal  regions  of  the  cylindrical  fibers. 
Combining  cohesive  failure  and  material  fracture  allowed  multiple  damage  mechanisms  to  be 
examined  at  once.  Crack  locations  in  the  combined  simulations  were  consistent  with  those  that 
only  included  matrix  fracture  without  cohesive  failure  behavior. 

One-ply  composite  cylinders  were  analyzed  as  a  means  to  isolate  the  influence  of  axial 
fibers  on  the  structural  capability  of  the  composite  cylinder  without  reinforcing  circumferential 
fibers.  For  the  one -ply  simulations,  fracture  occurred  at  significantly  lower  external  pressure  pe. 
Matrix  fracture  occurred  similarly  along  the  fibers  but  then  propagated  in  the  ^-direction. 
Considerable  damage  occurred  along  the  fiber/matrix  cohesive  interface  as  well. 

This  research  demonstrates  the  ability  to  model  multiple  modes  of  failure  in  composite 
materials  using  XFEM  and  cohesive  interface  failure — thus  improving  modeling  efficacy, 
diminishing  the  need  for  build-test-build  methodology,  and  vastly  reducing  cost  and  time 
requirements. 
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